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Abstract. We show that, if the long GRBs are produced by the collapse of massive 
stars, GRBs and their afterglows may provide a powerful probe of cosmology and the 
early universe. 



INTRODUCTION 

There is increasingly strong evidence that gamma-ray bursts (GRBs) are asso- 
ciated with star- forming galaxies [1,2,3,4] and occur near or in the star- forming 
regions of these galaxies [2,3,4,5,6]. These associations provide indirect evidence 
that at least the long GRBs detected by BeppoSAX are a result of the collapse of 
massive stars. The discovery of what appear to be supernova components in the 
afterglows of GRBs 970228 [7,8] and 980326 [9] provides tantalizing direct evidence 
that at least some GRBs are related to the deaths of massive stars, as predicted 
by the widely-discussed collapsar model of GRBs [10,11,12,13,14]. If GRBs are 
indeed related to the collapse of massive stars, one expects the GRB rate to be 
approximately proportional to the star- formation rate (SFR). 



DETECTABILITY OF GRBS AND THEIR 
AFTERGLOWS 

We have calculated the limiting redshifts detectable by BATSE and HETE-2, 
and by Swift, for the sixteen GRBs with well-established redshifts and published 
peak photon number fluxes. In doing so, we have used the peak photon number 
fluxes given in Table 1 of [15], taken a detection threshold of 0.2 ph s _1 for BATSE 
and HETE-2 and 0.04 ph s" 1 for Swift, and set H = 65 km s" 1 Mpc" 1 , Q m = 0.3, 
and Q\ = 0.7 (other cosmologies give similar results). Figure 1 displays the results. 
This figure shows that BATSE and HETE-2 would be able to detect half of these 




FIGURE 1. Cumulative distributions of the FIGURE 2. The best-fit spectral flux distri- 

limiting redshifts at which the 15 GRBs with bution of the early afterglow of GRB 000131, 

well-determined redshifts and published peak as observed one day after the burst, after 

photon number fluxes would be detectable by transforming it to various redshifts, and ex- 

BATSE and HETE-2, and by Swift. tinguishing it with a model of the Lya forest. 



GRBs out to a redshift z = 20 and 20% of them out to a redshift z = 50. Swift 
would be able to detect half of them out to redshifts z = 70, and 20% of them 
out to a redshift z = 200, although it is unlikely that GRBs occur at such extreme 
redshifts. Consequently, if GRBs occur at very high (z > 5) redshifts (VHRs), 
BATSE has probably already detected GRBs at these redshifts, and HETE-2 and 
Swift should detect them as well. 

The soft X-ray, optical and infrared afterglows of GRBs are also detectable out 
to VHRs. The effects of distance and redshift tend to reduce the spectral flux in 
GRB afterglows in a given frequency band, but time dilation tends to increase it 
at a fixed time of observation after the GRB, since afterglow intensities tend to 
decrease with time. These effects combine to produce little or no decrease in the 
spectral energy flux F v of GRB afterglows in a given frequency band and at a fixed 
time of observation after the GRB with increasing redshift: 
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where L u oc v a t b is the intrinsic spectral luminosity of the GRB afterglow, which we 
assume applies even at early times, and D(z) is the co moving distance to the burst. 
Many afterglows fade like b ~ —4/3, which implies that F v (u, t) oc D(z)~ 2 (l+z)~ 5 / 9 
in the simplest afterglow model, where a = 26/3 [16]. In addition, D(z) increases 
very slowly with redshift at redshifts greater than a few. Consequently, there is 
little or no decrease in the spectral flux of GRB afterglows with increasing redshift 
beyond z ~ 3. 



FIGURE 3. Cosmological context of 
VHR GRBs. Shown are the epochs of re- 
combination, first light, and re-ionization. 
Also shown are the ranges of redshifts 
corresponding to the "dark ages," and 
probed by QSOs and GRBs. 
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In fact, in the simplest afterglow model where a = 26/3, if the afterglow declines 
more rapidly than b ~ 1.7, the spectral flux actually increases as one moves the 
burst to higher redshifts! An example of this is the afterglow of GRB 000131. Its 
peak flux F pea k was in the top 5% of all BATSE bursts and the break energy -Ebreak 
in its spectrum was 164 keV, yet it occurred at a redshift z = 4.50. We have cal- 
culated the best-fit spectral flux distribution of the afterglow of GRB 000131 from 
[17], as observed three days after the burst, transformed to various redshifts. The 
transformation involves (1) dimming the afterglow, (2) redshifting its spectrum, 
(3) time dilating its light curve, and (4) extinguishing the spectrum using a model 
of the Lya forest (for details, see [15]). Finally, we have convolved the transformed 
spectra with a top hat smearing function of width Av = 0.2u. This models these 
spectra as they would be sampled photometrically, as opposed to spectroscopically; 
i.e., this transforms the model spectra into model spectral flux distributions. 

Figure 2 shows the resulting spectral flux distribution. The spectral flux distri- 
bution of the afterglow is cut off by the Lya forest at progressively lower frequencies 
as one moves out in redshift. Thus high redshift afterglows are characterized by 
an optical "dropout" [4], and VHR afterglows by a near infrared "dropout." We 
conclude that, if GRBs occur at very high redshifts, both they and their afterglows 
can be easily detected. 

GRBS AS A PROBE OF COSMOLOGY AND THE 

EARLY UNIVERSE 

Theoretical calculations show that the birth rate of Pop III stars produces a 
peak in the SFR in the universe at redshifts 16 ^ z % 20, while the birth rate 
of Pop II stars produces a much larger and broader peak at redshifts 2 < z ^ 10 
[18,19,20]. Therefore one expects GRBs to occur out to at least z pa 10 and possibly 
z ~ 15 — 20, redshifts that are far larger than those expected for the most distant 
quasars. 

Figure 3 places GRBs in a cosmological context. At recombination, which occurs 
at redshift z = 1100, the universe becomes transparent. The cosmic background 
radiation originates at this redshift. Shortly afterwards, the temperature of the 
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FIGURE 4. The cosmic SFR 
Rsf as a function of redshift z. 
The solid curve at z < 5 is the 
SFR derived by [25]; the solid 
curve at z > 5 is the SFR calcu- 
lated by [18] (the dip in this curve 
at z « 6 is an artifact of their nu- 
merical simulation). The dotted 
curve is the SFR derived by [24]. 



cosmic background radiation falls below 3000 K and the universe enters the "dark 
ages" during which there is no visible light in the universe. "First light," which 
occurs at z 20, corresponds to the epoch when the first stars form. Ultraviolet 
radiation from these first stars and/or from the first active galactic nuclei re-ionizes 
the universe. Afterward, the universe is transparent in the ultraviolet. 

QSOs are currently the most powerful probes of the high redshift universe. GRBs 
have several advantages relative to QSOs as probes of cosmology. First, GRBs are 
expected to occur out to z 20, whereas QSOs occur out to only z pa 5. Second, 
very high redshift GRB afterglows can be 100 - 1000 times brighter at early times 
than are high redshift QSOs. This makes possible very sensitive high dispersion 
spectroscopy of the metal absorption lines and the Lyman a forest in the spectrum 
of the afterglows. Third, no "proximity effect" on intergalactic distances scales is 
expected for GRBs and their afterglows, in contrast to QSOs. Thus GRBs may 
be relatively "clean" probes of the intergalactic medium, the Lyman a forest, and 
damped Lyman a clouds, even in the vicinity of the GRBs. 

The important cosmological questions that observations of GRBs and their af- 
terglows may be able to address include the following: 

• Information about the epoch of "first light" and the earliest generations of stars 
from merely the detection of GRBs at very high redshifts; 

• Information about the growth of metallicity in the universe in the star-forming 
entities in which the bursts occur, in damped Lyman a clouds, and in the Lyman 
a forest from observations of the metal absorption line systems in the spectra of 
their afterglows; 

• Information about the large-scale structure of the universe at VHRs from the 
clustering of the Lyman a forest lines and the metal absorption-line systems in the 
spectra of their afterglows; and 

• Information about the epoch of re-ionization from the depth of the Lyman a 
break in the spectra of their afterglows. 

Below we consider the first of these questions: the epoch of "first light" and the 
earliest generations of stars. 
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FIGURE 5. Left panel: The number AT* of stars expected as a function of rcdshift z (i.e., the 
SFR from Figure 4, weighted by the differential comoving volume, and time-dilated) assuming 
that CIm = 0.3 and f^A = 0.7. Right panel: The cumulative distribution of the number of 
stars expected as a function of redshift z. Note that w 40% of all stars have redshifts z > 5. The 
solid and dashed curves in both panels have the same meanings as in Figure 4. 

GRBS AS A PROBE OF STAR FORMATION 

Observational estimates [21,22,23,24] indicate that the SFR in the universe was 
about 15 times larger at a redshift z ~ 1 than it is today. The data at higher 
redshifts from the Hubble Deep Field (HDF) in the north suggests a peak in the 
SFR at z w 1 — 2 [24] , but the actual situation is highly uncertain. 

In Figure 4, we have plotted the SFR versus redshift from a phenomenological 
fit [25] to the SFR derived from submillimeter, infrared, and UV data at redshifts 
z < 5, and from a numerical simulation by [18] at redshifts z > 5. The simulations 
done by [18] indicate that the SFR increases with increasing redshift until z « 10, 
at which point it levels off. The smaller peak in the SFR at z 18 corresponds to 
the formation of Population III stars, brought on by cooling by molecular hydrogen. 
Since GRBs are detectable at these VHRs and their redshifts may be measurable 
from the absorption-line systems and the Lya break in the afterglows [4], if the GRB 
rate is proportional to the SFR, then GRBs could provide unique information about 
the star- formation history of the VHR universe. 

We have calculated the expected number N* of stars as a function of z assuming 
(1) that the GRB rate is proportional to the SFR 1 , and (2) that the SFR is that 
given in Figure 4 (see [15] for details). The left panel of Figure 5 shows our results 
for N*(z) for an assumed cosmology VL M = 0.3 and = 0.7 (other cosmologies give 
similar results). The solid curve corresponds to the star- formation rate in Figure 4; 

^ This may underestimate the GRB rate at VHRs since it is generally thought that the initial 
mass function will be tilted toward a greater fraction of massive stars at VHRs because of less 
efficient cooling due to the lower metallicity of the universe at these early times. 



the dashed curve corresponds to the star- formation rate derived by [24]. Figure 5 
shows that N*(z) peaks sharply at z 2 and then drops off fairly rapidly at higher 
z, with a tail that extends out to z ~ 12. The rapid rise in N*(z) out to z ~ 2 
is due to the rapidly increasing volume of space. The rapid decline beyond z ~ 2 
is due almost completely to the "edge" in the spatial distribution produced by the 
cosmology. In essence, the sharp peak in N*(z) at z « 2 reflects the fact that the 
SFR we have taken is fairly broad in z, and consequently, the behavior of N*(z) is 
dominated by the behavior of the co-moving volume dV(z)/dz; i.e., the shape of 
N*(z) is due almost entirely to cosmology. The right panel in Figure 5 shows the 
cumulative distribution iV*(> z) of the number of stars expected as a function of 
redshift z. The solid and dashed curves have the same meaning as in the upper 
panel. Figure 5 shows that for the particular SFR we have assumed, ~ 40% of all 
stars (and therefore of all GRBs) have redshifts z > 5. 
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